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INFORMATION THEORY AND THE EARTH’S 
DENSITY DISTRIBUTION 

David Parry Rubincam 
NAS-NRC Resident Rcsearcli Associate 

ABSTRACT 

The present paper argues for using the information theory approach 
of Jaynes (1 957) as an inference technique in solid earth geophysics. 
A spherically symmetric density distribution is derived as an exam- 
ple of the method. /\ simple model of the earth plus knowledge of 
its mass and moment of inertia leads to a density distribution which 
is surprisingly close to the optimum distribution of Bullen (1975), 
Future directions for the information theory approach in solid 
earth geophysics as well as its strengths and weaknesses are discussed. 
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INI’ORM Allow nil ORY AND Till I AlUll’S 
Dl NSri Y DISTRIIU'TION 


1.0 INlRODl'niON 


\Vc wisli tt- introiUico to solitl eaith iicopli> sics a mothoil ot scicutil'ic iniorona* which 


has had jircat success in statistical mechanics (see, e.^i,, lay nest 0^57, liiluis(l*>nl); 


Kat/ (l%7); and Baierlein (l‘)71 D and in spectral analysis (e.g. Burg Smylie iM 

ah (l‘)7.D; and Oraber (O)7(0>. It is the inloimation theory approach of Jaynes il‘).S7). 
based upon Sliann^in’s ( 1 ‘MS) infonnation measure. We will illustrate the approach by 
inferring a density listribntion for the earth based on knowledge of its mass atul moment 
of inertia. The earth is assumed to be splieiieal and the density distribution spherically 


symmetric. 

riie nature of the inference problem is the following. We desire to know '’,hat the 
density distribution p(r) is as a function of radial distance r from the center of tlie eartl',. 
Suppose the only information we have is us mass Mj,- and moment of inertia both ol 
which depend upon jiKr). Clearly we do not have enough information to say what the 
density distribution pfr) actually is. Any proposed distribution whicli satislies tlie mass 
and moment of inertia is uommique; there are infinitely many other distributions which 
also satisfy the given data. 

I'here are several methods for dealing witli this problem, (l or a general discussion see 
Bullen (l‘>75, pp. (>0-(>4).) The approach of Backus and (lilbert (l‘>(i7, l')o8) is to study all 
solutions consistent with the given data; this is called the geophysical inverse problem. The 
Backus-Gilbert approach has been used extensively. See, for example, Gilbert et al . (1^7.)); 
Parker (I ‘)77a, l‘)77b); Jordan and Trauklin (l‘)7n;and references cited by Parker (1‘) 77 a. 
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I*>77h), Richiinls ( l‘)75), Andoison (l‘)75), iiiul l-ngiliilil el ;t|. 1 1‘>75). Press ( 
iiilopted a Monie (‘arlo leelinique of lesiinji a wide range of models against the data and 
retaining only those whieh agreed with it. However, ilte eommonest metlunl by far is model- 
ing: by introducing other assumptions, the answer Incomes unique, I'lie assumption of the 
Adams-Williamson equation, for instance, plus the known mass, moment of inertia, seismic 
velocities, and surface demsity determine a unuiue density distribution {Alierman el al. 1 
pp. SO-Sl). Of course a difficulty with this approach is that the assumed conditions may 
not hold. 

Suppose we look ; l the (iroblem fnun the following viewpinnl. If we had to pick one 
answer (in onr case, density distribution) from all the possible answers which fit the data, 
which would we pick as the most likely? I'o put it bluntly . what is our "best guess"? It is 
of extreme interest that the information theory approach of Jayues(l‘)57> provides an 
answer to this question. (Haiericin 1 1 ‘>71 ) has an excellent general discussion of the infor- 
mation theory approach.) 


2.0 INl-ORMA PltW Ml A.SURb 

At the heart of the approach is .Shannon’s ( I nqs) information measure 


Ml (I’l . IN 


..Pn) 




PdnPi 
i-1 ‘ ' 


( 2 . 1 ) 


Here Pj is the probability that the ith of N possible answers is true and K is a positive con- 
stant. riiis function was originally termed the entropy function ( Tribus and Meirvine, I ‘>71 , 
p. I SO), due to Us similarity to thermodynamic entropy. I'or this reason the information 
theory approach is often caHai the maximum entropy method, or MPM for short, riie rela- 
tionship between the information measure and thermodynamic entaipy is deep, but the 


l\\‘0 .lu' Ui't iili'iUu'.ii ( M-iu'i Iv'iii, ‘I . I'p. ’I *,t ■} ’SV !«' I'lX.uiiMi'ii w i‘ Will 

liau'iK'm p (>-1' .iiui ».ill S!i.tiiiu’n\ jiiioi m.iliou iiummiu’ Mld’i . T*. . 1\1. 

wlu'io Ml si.uiiK l«M "Mi'.miii' liiii'im.iti.'u", Hi, .mu'iuti iM iiiii'inuthm iK'i\ii\l Jo ilo 
toimuu' whu'li answiM is vOiuM 

\VV Will ih'J piiWv* tll.lJ \|| IS ,i MtilsuK' ol MIIsmiu; iiiI,m iiulh'ii . piOv'Is .lie p.iNi'tl l'\ 

Sh.timon I l'>"l.s'> ,nul U.iieilein 1 1‘* 'll |\.(j|iot, we will metel\ iiulu.iJe its pl,iiisilMlit\ with 

,iu e\.imple. Mill litsj we note tli.it Ml • 0, the .nnouut ot mloiui.itioii needoil to siiijtle 

out the ks'iu\'t .luswei IS novel noiMtive 1 his is i,'ii.iinl\ .in intmtivelv ilesu.iHe pu'peitv. 

Now let us suppose ih.ii .ill lU the I'loivihiliiies .ue .■.pi.tl, lu this ease Ml attains its nia\ 

iiluiiu value. 1 Ills avw'iils with tniinlton we aie siiielv in ,i stale ,'l 'ii.ixununi ifinotanee 

li.e, luvil the most uilv'ini.ilioiH it we ean t.ivoi no .iiiswei ihove aiU'lliei in leiins of 

pi oKihilttv . Suppose n,'w v\,' h.ive ,iisev'veie,l lh.it die ith i'ossiImIiIv is the ev'iieel answei. 

I hen I’l I ,iiul I', v' toi i ^ t !l,»w tmt.Ii intoim.iiion is missuu! now ' In this ease 

r,Inr, Hill O, ami r.lnl* 0 loi I 1 dw viuue of lim \!n v til. Ihus Ml 0 no 

\ *0 

inloiination is nussiue, we h.ive the .iiiswei, I his .ilso .u.ouls with intuition. Nv'iuialh 
out tpou'ianee lies I'eiweeii tliese two eviiemes. .uul Ml takes on values aeaniluiiilv he 
tvveen its iii.iMinuin .uul th lienee Ml is .i pl.iusihle uie.isuie I'l nussuii’ luroinuitioii. 

We slioulil point out heie th.il Ml is not \hniensioiiless il iliiiuiulson. piiv.tte eoinnui 
uieaiiou, 1'1‘uKa taet tlui iloes not .ippeii to tv evplieith noteil in luhusilO(,n, 
ll‘>(vH, oi Haieilem iH’HV It eaities units ot intoiiiutu'n. I oi example, it we ehan>*e 
the h.ise ol the to>'.,iiilhm lU e>|. t '.n tiom e to wliieli eli.iue.es K. tv' .1 new ev’iist.int K', 
atui set K. 1, then Ml l'l’,lot^>P, .uul Ml is measuieil m hps. In the lotlowinn ileveK'p- 
ment we will lel.iin (he uatmal loe..uithm base aiul set k I . so that Ml is measmeil ir, 


iiats urom natural units). Wo will supjuoss tho units m tho tollowinji ilovolopniont. but 
it slunilil bo lomomboivil that Ml is not a iliinotisionloss quantity. 

JAVNt' .S* MKlNr i n I t)l' MINIMI’M IMU Jl'niH' 

riio ossoiKV ot tho inlonn.ition thovuv appiv'.ioh is this, ohooso tho pioKibilitios P| , I’l. 

. . . ol tlu- posNiblo outoonios to inako Ml as latno as possiblo. subioot to tho oonstiaints ol 
tho known data, lliis is J.iynos* piinoiplo ot nuuimum ptojiulioo iliibus aiul Rossi, 
llonoo tho intormatuMi thooiv appio.ioh is a rational inothod tor assii’ninn piobabililios. 

lot us illusliato tho toohiUi|Uo with an osainplo. .Supposi,* that wo do not know tho 
mass ot tho oailh oxaoih . but iduo to oxpoiimontal ortoi, sa\) it must bo ohoson iVom 
tho valuos Mj. M>, . . . M^. \sido ftom 1. this is all wo know. Wo must tind Tj. 

tho piobability that Mj is tho ooiiooi mass, by inaMini.'inj’ Ml, I’liis is dono by lakinn 
tho partial donvalivo ol 

With rospoot to oaoh l*j and sottiU}’ it oqual to .'oio. Iho o^^ is a I aniaiigo muliiplioi whioh 
iiisuros that all ot tho piobabililios .uid ui' to 1. (.’ariyinj’ out tho pioooss yiolds 

-InP, - I t ivq 0 
or 

\\ constant 

Iho unknown op nuy bo lonnd tiom tho constraint 

N 

i: I’i 1 

i 1 ' 

giving 


■I 




i\ -- I N 

All of till’ prob.ibililios .iri' oiiiul, Wc >^nv'\v nothini! about tlu* vatjous Mj to favor osk' 
particular value over another. 

Now su)^pv^sl• we obt.riu further mioriuatiou; we learn Miat the espectation value 
of the mass is 

N „ 

i: RM, Mp 
it ' ‘ 

We reassign probabilities in accordance with Jaynes’ principle. 

•M-nViil’i + oo“l'i + iMi:i\Mj) l> i - 1. ’ N 

dl’i 

giving 

i 

where iV 0 and 0| are lagr.urge luultipheis to be found from the eimstraints 

I’P, 1. I'l’jMj Ftp , 

Note that om method is comt'letelv analogous to tliat of the canonical ensemble in sta- 
tistical mechanics. Indeed, tire irra are nra tics ate identical. The only difference is in lire 
philosophical basis, which indicates that the method has broad .ipplicabilit^ and is not 
confined *mly to statistical ir’’chanics. 

Obvionslv assigning probabilities is rather easy. Hut clearly lire probabilrties do irot 
roprosenl frcduencies; their vairres change the mourent we aciiuirc new iirtornration. What 
is the point in ma\inti/ing Ml ’ Why not do sotrrething else'.’ 

lire users of Jaynes’ prirrciple have a povverfn! ,ngnnrent in its favor. I'lie probabil- 
ities are iirdeed not nrcre frciiuencios, they say, A probability rcpreseirts the “degree of 
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rational bi'liot” (Haioilcin. l‘>7l, p. l."i tliai a paitkular answer is correct, a more general 
notion than a IreMiiency. l5i,e l’o\ (hMo. for the quantitative basis for this view 
of probabililv .) And since Ml measures the anonnt of information needed to determine 
the correct answer, miy melhoil for iissigiiing probiiblHties which does not mn\inii/c Ml 
imdei* known eons' i. lints (knowledge) taeitly assumes information it hasn’t got! In other 
words, if someone assigns probabilities not in acandance witli Jaynes’ principle, that 


ivrson is prejiidi. mg the probabilities without foundation in the known data. I’hus the 
name, "principle of minimum prejudice." 

riiis point is particularly cleat in the evami'le where we knew one of the Mj was 
the correct annwer. but h.ul no other information (other than i-'l*, 1). In this case 

Jaynes' principle assigns equal probabilities to all outcomes. We aie completely ignorant 
as to which itnswer is correct. If someone uses some otlier principle, ,md assigns (say) a 
larger probability to M| than to the other Mj, we can say, "^ou taviucil Mj as being 
the most likely mass over all of the others. What b.isis (i.e. information) do you have to 


ilo that”’’ While the argument is powerful, the information theory approach is not w'itli- 
out its problems, Wa' will discus> some ot these later. 


4.0 INFORMATION 1111 OKV orNSU Y ntsrmtu'TioN 

W'e are almost in a position to find the information theory density distribution in- 
side the earth, using knowledge of the expectation values of the mass and moment of 
inertia. The only thing left to do is set up the problem. W^* whll make heavy use of the 
methods of statistical mechanics; p.irticularly, the grand canonical ensemble. 

Imagine a three-vlimentional ('artesian coordinate system with its origin at the center 
of the earth. The grid sy stem will divide up the earth into many cubes of identical 


0 



\ oliiino V, just .is ijmi'Ii i'.i|vi di\ iilos up ,i pLiiu' iiUo viu.iri‘>. i-qiLil .ium. We 


i\m .ipproxim.iti’ tlu* spluMu'.il sun'.uv oi the e.itth ,is eli’seh .is we like l'\ m.tkm^ the 
eiihes ,is sni.lll .is we like, let he the \eetoi lu'in the eeiltei ol the e.irtli to the itli 
eiihe .mJ set I r[ I i(. I et the iii.iss oi the e.nth he the sum ol the iiusses ol' a l.ii^te 
miiuhei ol iiulistinsuisluhle p.utieles. e.ieli w ith m.iss m. I et there he ii, p.iitieles m tiie 
ith eiihe. The m.iss M|, aiul moment ot inertia (’| ol tlie earth are then 


Ml, i;n,m 

I 

Cl : ,?i:n,inrj* 
1 


i4.n 


where the suhseript i runs over .ill tlie etihes eompiisiiif> the earth. The laelor t ’ . 1 ) in 
the seeotul equation m ikes use ot tlie *'aet that the density distrihution is spherieallv sym- 
metrie, .nut takes eare of ij hemu the dislanee Itoin the eentei ol the earth I'l a euhe and 
not the dist.inee to some axis of loiation. 

1 et us lem.irk heie that we ha\e ehosen euhes ol ei|u.il vi'lume so as to treat all 
rejjioiis ot the e.nth identieally, and indislinputsli.ihle p.utieles heeanse the interehanjjini; 
ol p.irlkles leaves the densitv ilislrihulion unalleeted. We make no eomniitment as to the 
values of m and ,\s we shall see. ihev drop out of the final eiiu.ition for the density 
disfrihution. 


We .ire ready to hepin. Our information will he that the earth may he made up of 
any numher of ixirtioles, hut that the expeol.ition value of the mass l'l\Mj and moment 
of inertia l'l’,Oi are known to he M|. and r*|.'. lespeetivelv . In praetiee, Mp and Op will 
he the experimentally determined values. W hat we will do is find hj, the ex]>eetation 
value for the numher of p.utieles in the ith euhe. I'he prohahitities are eomptited aceord- 
inp to j.i\ lies’ prineiple ol minimum preiudiee; 


*T 










wlunc 


/ v^.‘'iM| + *ii ’ 1 1 
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V'rom c\i, I’l.n we nuiN wute 


M, l‘n| 


tn 


2m I 




i 4 ..n 


where tire subseni't i on the n, h.»s been snpinesseil. Vlre problem irow looks e\.iell\ like 
that ot the jnaml e.inoiueal ensemble, with n, plaving the role of oeenpation n'linbeis, tj* 
the role v’f enetps lewis, .m\l e^i. i4.2l the piaiul p.utitron fnnetion. The treatment of 
this problem sna\ be found m an\ statularil statistieal meehanies text. We elux'se to lol 
low Morse tl'ho), pp, .t22-d2o\ 


Vsing evt. pl,.M in eo. t>t.2V we have 




t4.4\ 


where we have ledefmeil oj in as oj and 2 di\_jni aso^'. Note that 

din/ 
do I 


r4,sf 


.otin. + ajiiiji,- 

^ .TV 




a result that we will make use of shortlv . 


i #»*'• f 
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Let us now rewrite eq, (4.4) ;is a summation over the possible values ol Hj insteail of 
over i, For elistiiitfuisliable particles it is (Morse, h)o‘). p, ,)24) 

-y - V Ni +aii.n,r,“ 

“ I, , 1 ,1-1 ^ 1 l 

with N = H| + ii 2 + ... . 

The thint; to do now is make 7 . matiiematically tractable. We do this in the follow- 
ing manner, Assume the cubes are so small that the chances of two particles sharing the 
same cube are negligible. This is equivalent to assuming the particles follow Maxwell- 
Boltzmann statistics. To take care of the indistinguishability tif the particles we can then 
divide the above etpiation by N! ami obtain 

111,112... ntinal 

as the ; r,' '’ mate value for the grand partition fu.r.iion. I'lirther. since 

I 


I 3E 


we can separate Z into factors for each cube; 




where 


by virtue of 


7 = v,.(ftt +a2rj‘)ni , v,,(oi + ft^r ‘)n2 
«1 «2 

= Z| ‘Zs* . . , 


Zj = exp(e“J 


e’^ = S x'Vn! 
n=0 


From cq. (4.5), and above we have 

dInZ 


3a 


= ve“'e“2^r = £n, 

1 j j J 
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whore evidently 



The density distribution is obviously 

By the assumption of sphericul symmetry tor the density distribution we can drop the 
subscript and write 

p^(r) « p(0)e“2^^ (4.6) 

which we will take ns the desired information tlieory density distribution. The two constants 
p(0) =-^e“l and may be found from our knowledge of the expectation value for the 
mass and moment of inertia: 

Me = 4^*** p(r)r2dr = 5.976 X gin (4,7) 

Cjj ** ~y^***^ p{r)r^dr = 8.068 x 10^*^ gmem^ (4.8) 

where be is the radius of the earth and our numerical values have come from Stacey 
(I960, p, 279). hi eq.s (4.7) and (4.8) we have assumed that the cubes are so small that 
we may switch from summations to integrals without serious error. By numerical inte- 
gration of eq.‘ (4.7) and (4.8), we find tliat 

p(r) = 12.30c''* gm/cm^ (4.9) 

is our “best guess*’ for the density distribution based on the given data. 

A plot of eq. (4,9) appears in Fig. 1 , along with the “optimuni” density distribution 
given by Bullen (1975, p. 361), wiiich presumably gives the most plausible distribution on 
the basia of all the known data. Tlic two curves agree remarkably well, in view of the 
fact that the information theory density distribution makes use of only two basic pieces 
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ot iliUa lUiiss aiul moment ot mevtia. No seismie ‘.latii or eiiuations oi slate have Iveen 
mohuled in om inloinuUion. 

5.0 1 iniMU- IMKHTU^NS 1 Ml- rUl ORV 

Vho lesnll obtained above may be e.isity genernli/ed to melnde .my known volume 
integrals ol llie vlensity ilisU Unit ion. Supposing that there are I sneli integrals having the 
Imin ^ 

J f>tr*)l ,u^ilv ■ I j U 1. ... I > 

wliime 
of eiiith 

the resulting average density distribution is 

; oir^l eonsl • espiojliu*^^ * • . . . 

I'he I agrange multipliei^ Oj are tound tiom the known values I'j. Note that the alnne 
lesiilt IS not lestiieled to the spherually symmolrie ease. Uesides the mass and moment 
of inertia, the splieiie.il harinonie eoellieients and Sj„^ of the earth’s gravitational 
potential immediatelv eonie to mind as integrals having this form. We mleiul to publish 
lire lesultmg pii*^l based or, the gravity iietd eoellieients in the near 'oture. 

riie ne\t obvious extension of the theory is to assume the e.iith is an elastie body 
so as to ineliide the ekistie pammeleis fui^l and Mr^l in addition to the ilensitv distubu 
lion ptr*^> .1'' unknown nuantUies to estnuale. This will allow seismie tiavel limes, free 
oseill.ition periods, and body tide obseiv.itions to used, all of whieh depend upon fifi*^!, 

\(i^f and (Jrabei 1 1 ‘> ’ n has made a start in this direetion using torsional vibrational 
modes of the earth, but mueh work needs it> be done. Inelusion ot ehemieal eoinposition 
atiil equations of state .ire lUher j'ossible direct kms for the theory. 

The goat here is to put into the pioMem all of the physics and data that we know 
and maximi/e the remaining missing information. If we s’oine to a point where the phvsies 
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iiml data iletcrminc a unique answer, there is no need lor information theory or any otlier 
inferenee teeimiq\ie. There is no remaining uncertainty; we have the answer. Since this 
liappy state of affairs is unlikely to occur soon, geopliysics has a real need for sound 
methods of infereiice, 

O.O OISCUSSION 

We will discuss hero ja»me problems and questions about the information theory ap- 
proach. The two main problems which face it have been discussed extensively in the 
literature and only bear indirectly on the geophysical problem at hand, Hence we will 
merely point them out aiul then move on to questions more relevant to geophysics. 

The first problem is that the information theory approach apjwars unable to deal with cer- 
tain kinds of information. For example, suppose we tlipa coin 100 times and it comes up heads 
75 times, (’learly we have information on whether the next fliji will be u heads or a tails. But 
how do we maximi/e eq.i2.h using this data? Information theory seems to be silent on this 
question. This coin flip problem is discussed by Rowlinson ( l‘>70) and references he cites. 

Tlie second question is what to do about continuous probability distributions. It has 
vexed even the most ardent proponents of the information theory approach. 

Shannon (, l‘^4S. p. 028) proposed as the appropriate genera li/ation for the continuous 
case in one dimension the function 

Ml - “K / ptxllnptx'ldx 

••■OO 

where p(x) is the probability distribution and \ is a continuous parameter. One difficulty 
is immeiiiately apparent. It takes the logarithm of a dimensioneil quantity: p(x> has 
dimensions of inverse length if x is a length. Siiy. The biggest stumbling block, however, 
to the use of the above equation is in assigning prior probabilities. To ilhistrate. suppose 


X roprcsi’uls tlio spcetl of a partiolt' amJ is known lo lie between values Xj and Xi. We 
fiave no further information. If we apply Jaynes* prineiple we find p{x) * constant for 
X( < X and zero outside tlje interval. But if we had taken the kinetic energy of the 
particle {which depends on x~) as the continuous parameter and found the probability 
distribution as a function of kinetic energy, we would have again found it to be constant 
over the interval. The two are inconsistent: a constant distribution for speed implies a 
nonconstant distribution for energy and vice versa. Thus the problem, Jaynes (l‘J63, 
l‘)08) and Holxsoti and Cheng (I‘f73) argue that the above equation must be modified, 
while Tribus and Rossi (I ‘>73) and Batty (1 ‘>74) feel that it is the correct equation and 
that the inconsistency is a pscudoproblem. Rowlinson (l‘)70) also discusses this matter. 

Getting back to the ilensity distribution derived here, a question arises; why assume 
that the cubes are sparsely occupied, thus giving Maxwell-Boltzmann statistics? This 
would seem all right for a dilute system, but in the case of the earth it would appear more 
reasonable not to limit the number of particles in each cube. 

If in fact we allow an unlimited number of particles to occupy each cube, then we 
are dealing with Bose-Hinstein statistics, I’inding Uj in the usual statistical mechanical 
fashion (Morse, 19o‘>, p, 32b) yields 

1 

h(rj) = 1 

J g«l + «2fj‘ _ 1 

using the mass and moment of inertia as the eon.‘:t»'a fists so that 

■f 


^r) = 


e«l + *2r‘ _ , 


(0.1 J 


Here we have a problem. There arc more unknowns than constraints. If we knew how to 
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choose -1^, then we could litui iV| .md a 2 Uie constraints of muss ami moment of 
inertia. Unfortunately, we have no clear guidance in this matter. 

rhcrc is a way. however, to neatly sidestep the problem. We introduce a third piece 
of information: we assume we know /flap), the v due oi the density at the earth's sur- 
face. Using this information in cq. (tnl) yields 

and we use our knowledge of Mj. and U to find the two multipliers, faking ^tapO as '.S4 
gm/etU'^ to agree with the surface density of rocks iStacey. p, 104), we timl by 

nuincric-al integration 


/•5U) 


P(i) , K,;/, ~ I"” "'”' 

0.44e‘’**‘‘ *'t + 0.0007 


U'.J) 


This distribution is plotli-d in Tig, ^ along with llulleiTs tl07.s> optimum distribution. 
1‘his Hose-binstein distribution is almost indistinguishable from the MaxwelbBolt/rnann 
distribution of rigure I amt the distinction between using the two is academic in this 
instance, at least. 

So we could use llose-l'instcin statistics in the information theory approach to the 
earth's density distribution. While it is probably conceptually superior to Ma\welM)o!t/- 
tnann statistics, it is also a little more eomplicated mathematically, as a comparison of 
eq,s 10,1) and (O.') shows. The Ma\welMtoll/mann case in future investigations should 
probably be investigated first, being simpler. 

Anotiier question arises about the information theory approach. Suppose we are 
again discussing the example of the earth’s mass, where we have N masses to choose from 
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anil know that i-'PjMt - Mu' inlovnwtion theory probability vbstributjon is c\iHOo“ • * 
i\(M, ^ as toutul earlier. Yet the esperimenia! probability vlistribution is probably a tJansMan 
with the form eonst • e\pt-(Mi-Mi.‘>“ 'o*> rather than the simple exponential Uistnbulion oJ 
mil - mation iheorx . with In'tli Mj.- ami n bemp known. The same may be s;hil for the mo- 
ment of inertia. How ilo we haiulle this? 

It is aetiully rather easy to obtain a thuissian Uistribution from the information theorv 
approach. Assuming that we know 

we ilifferentiale 

-I'lVnrj + Ooi.Tj -oj 
with re.spect to e.ich I’j, pivinp 

p. . 

which is tlic vlcsircil tlaussian. Presumably limlinp nu]) iisinp a tlaiissian is ttiffic.iU, but the 
iiuestion is piob.il'iy moot, As tonp .ts we ileal with a narrow ranpe of values ot mass near 
Mp-.U ivouKl seem to matter little whether we use exponential or tians.sian ihstributiivns. 

As we liave seen, tlie inlormation theory api'roaeh lioes have its problems; but the 
stienpth.' of the approaeii are main . Its philosophieal basis is satisfy uip: there is no umvar- 
rantoil weighlinti of possible answers. It is rational auii obieetivo; every one usinp the approaeii 
will obtain the same answ er.s, given the s.ime il.ita. It gives the “best" answer on the b.isis 
of very lillte vlala. Hie mathematies of the theory is stamlaril" that of statistical meohanies. 
Observational errors have a natural I'laee in the approaeii. It provides an alternative to ex- 
tensive modeling, riiullv, its generality gives it great power; it ean be applieil eiiually well 
to tbe statistics of a pas or the interior of the earth. The information theory approaeii 
should find broad use in solid e.irtli geopliysies. 
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r(km) 


Figure \, The infonnatjoji theory ileusity distribution using Maxwell-Uolt/iuami statistics 
(curve A) anil the optinnun density distribution of Bullen ( I d75) (curve B) are shown as a 

function of radial distance r. 
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